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ABSTRACT 

New 1.4 GHz VLA observations of the pulsar-powered supernova remnant 3C 58 have resulted in 
the highest-quality radio images of this object to date. The images show filamentary structure over 
the body of the nebula. The present observations were combined with earlier ones from 1984 and 1991 
to investigate the variability of the radio emission on a variety of time-scales. No significant changes 
are seen over a 110 day interval. In particular, the upper limit on the apparent projected velocity of 
the wisp is 0.05 c. The expansion rate of the radio nebula was determined between 1984 and 2004, 
and is 0.014 ± 0.003% yr" 1 , corresponding to a velocity of 630 ± 70 km s _1 along the major axis. If 
3C 58 is the remnant of SN 1181, it must have been strongly decelerated, which is unlikely given the 
absence of emission from the supernova shell. Alternatively, the low expansion speed and a number 
of other arguments suggest that 3C 58 may be several thousand years old and not be the remnant of 
SN 1181. 

Subject headings: ISM: individual (3C 58) — supernova remnants 



1. INTRODUCTION 

The supernova remnant 3C 58 (G130. 7+3.1) was clas- 
sified as a pulsar wind nebula (PWN, also plerion; e.g., 
iWeiler fc Panagiall978|) long before a pulsar was known. 
Its pulsar, PSR J0 205+6449, was only recently detected, 
first in the X-rav |Murray et al.ll2002|) and then in the 
radio l)Camilo et al.l I2QQ 2|) . 3C 58 i s at a distance, 
D, of ~3.2 kpc l|Roberts et alJ I1993D 1 . 3C 58 is a 
bright radio source, and previous radio imaging observa- 
tions of the synchrotron nebula have shown it to have 
a center-brightened morphology, and a siz e of 6' x 9' 
llWeiler fc Seielstadlll97lt iWilson fe Weilerlll976T: iGreenl 
ll986URevnolds fc Allerlll988t iBietenholz et al.ll2001b[) . 

The pulsar in the center of the nebula is loosing rota- 
tional energy due to magnetic braking, at a rate which 
can be calculated from the pulsar's spindown rate, and 
is 2.7 x 10 37 erg s _1 ( assuming a neutron star moment of 
inerti a of 10 45 g cm 2 : lMurrav et aDl2002t |Ransom et al.l 
2004) . This is the third highest spindown power of any 
Galactic pulsar, the highest being that of PSR B0531+21 
in the Crab Nebula. This energy flows from the pulsar 
as a wind consisting of relativistic particles and magnetic 
field, which then powers the visible synchrotron nebula. 
The details of this energy flow, however, are not yet well 
understood. In the Crab Nebula, which is the best stud- 
ied pulsar wind nebula, rapid variability has been seen 
throughout the region near the pulsar at radio, optical, 
and X-ray wavelengths. In particular, a series of features 
called wisps propagate outward at projected speeds of 
orde r 0.3 c. The Crab nebula wisp s are visible in the ra- 
dio llBietenholz et alJfeOOlal 120041). i n the optical (e.g., 
iHester et allll996UTanvir et al.lll997l) . and in the X-ray 
( Mor i et al.ll2002j) 

These wisps are thought to be associated with the ter- 
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that the pulsar dispersion measure is a pproximately twice what is 
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mination shock in the pulsar's outflow, and indeed con- 
stitute almost the only observational window on that out- 
flow before it enters the body of the nebula. Radio fea- 
tures apparently associated with the pulsar outflow are 
visible in pul sar nebulae other than the Crab, for exam - 
ple in Vela llBietenholz et alJll991at iBock et aTl f2002). 
In particular jFrail fc MoffettJ l)1993|) showed that 3C 58 
also has a radio feature with a wisp-like morphology. It 
would clearly be of interest to determine whether 3C 58's 
wisp moves as rapidly as the wisps in the Crab, and to 
determine whether rapidly moving wisps are a general 
feature of young pulsar wind nebula, or a peculiarity of 
the Crab. 

Another outstanding question concerning 3C 58 is its 
expansion rate, in particular that of the synchrotron bub- 
ble. 3C 58 has often been associated with the super- 
nova of 1181 A.D. (SN 11 81: IClark fc Stephensor]ll977l 
iStephenson fc Greenlll999j) . implying an age of ~ 820 yr 
for the nebula. Such a young age suggests that the neb- 
ula should be expanding rapidly. The measured expan- 
sion speeds, however, o f both the synchrotron bubble 
([Bj£j£nh^lz_el al] |2QillbJ) and of the thermal filaments 
l|Fesenl Il983t iFesen et ail Il988l Ivan den Bergbl [l990), 
seem to be considerably lower than expected for an age 
of - 820 yr. 

Such slow current expansion speeds might be observed 
if 3C 58's expansion has been greatly decelerated, but 
this is unlikely for the following reason. The synchrotron 
nebula is expected to be expanding into the super- 
nova ejecta, which are themselves expanding. The syn- 
chrotron nebula is at fairly high pressure compared to 
that expected in the supernova ejecta, so in order for 
synchrotron bubble to have decelerated, the ejecta would 
also have to be decelerated. If the ejecta had been de- 
celerated, however, one would expect bright radio emis- 
sion from their interaction with the surrounding mate- 
rial. Such emission has not been observed despite nu- 
merous searches. (For more details on this argument, 
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see iBietenholz et alJl2001b|) . In addition, such decelera- 
tion of the ejecta would likely produce consid erably more 
therm al X-ray emission than is in fact seen ane et all 
2004). The low measured expansion velocities suggest 
that 3C 58 might be several thousand years old, and 
therefore not be the remnant of SN 1181. An accurate 
determination of the expansion speed of the synchrotron 
bubble is therefore important to reconstructing the his- 
tory of 3C 58, and determining whether it is in fact as- 
sociated with SN 1181. 

I obtained new radio observations of 3C 58 using the 
NRAO 2 Very Large Array (VLA), with the goals of ad- 
dressing both the question of possible rapid motions near 
the pulsar, and in order to obtain an accurate measure- 
ment of the present expansion speed of the synchrotron 
nebula. I describe the observations in § [3 discuss the 
results in § |3 and discuss my findings in § 0] 

2. OBSERVATIONS AND DATA REDUCTION 

I observed 3C 58 using the A, B, and C configura- 
tions of the VLA at 1.4 GHz. Details of the observing 
runs are given in Tabled Two separate A-configuration 
runs, separated by 110 d, were obtained in order to 
check for rapid changes at the highest resolution. To 
avoid bandwidth smearing, the observations were car- 
ried out in spectral line mode, with a total bandwidth 
of 44 MHz comprising 7 spectral channels, each of width 
6.25 MHz. The data were calibrated and imaged using 
NRAO's AIPS software package. The flux density scale 
and instrumental bandpass were determined using obser- 
vations of 3C 286. 

1 also re-edited the data of iRevnolds fc Allerl (j.988), 
taken in 1984 and also using the A, B, and C array 
configurations at 1.4 GHz, but without using the spec- 
tral line mode and having a bandwidth of 12.5 MHz. 
These data were then imaged and de-convolved in a sim- 
ilar manner as those from 2003/2004. Finally, I also re- 
reduced archival A-configuration 5 GHz VLA observa- 
tions of 3C 58 taken on 1991 June 30, w ith a bandwidth 
of 50 MHz (original results published by Frai l fc Moffetti 
1993). 

Obtaining high-dynamic-range images of 3C 58 
presents a challenge because the field around 3C 58 con- 
tains a number of background sources, including a nearby 
bright extra-galactic double. To obtain the best image of 
3C 58, these background sources must be included in the 
deconvolution, otherwise their sidelobes will degrade the 
image of 3C 58. Although maximum entropy deconvolu- 
tion is well suited for extended emission like 3C 58, the 
program available in AIPS (task VTESS) is not capable 
of treating the multiple fields necessary to also decon- 
volve the background sources. Furthermore, the decon- 
volution of compact objects like the background sources 
is better carried out using CLEAN rather than maximum 
entropy deconvolution. I therefore used a hybrid decon- 
volution scheme which made use of both maximum en- 
tropy and CLEAN, with the final images made as follows. 
First, a CLEAN deconvolution was performed on both 
the 3C 58 field and of six secondary fields to include the 
nearby background sources. Such multi-field CLEAN im- 

2 The National Radio Astronomy Observatory, NRAO, is oper- 
ated under license by Associated Universities, Inc., under cooper- 
ative agreement with National Science Foundation. 



ages were used to fully self-calibrate the data in both am- 
plitude and phase. For the final images of 3C 58, I sub- 
tracted the CLEAN components for all the background 
sources, in other words all the clean components except 
those for 3C 58, from the fully self-calibrated u-v data. 
The resulting u-v data set was then imaged and decon- 
volved using VTESS. The VTESS (maximum entropy) 
deconvolution reduces the rippling at small spatial scales 
which was visible in the CLEAN images of 3C 58 and is 
a known instability of the CLEAN algorithm. 

In order to determine whether there is any rapid vari- 
ability in 3C 58's radio emission, I made separate images 
from my 2003 July 7 and 2003 August 9 data. Time 
resolved radio observations of 3C 58 with the VLA, how- 
ever, present the problem that it exhibits structure at 
spatial scales from ~T0' down to <1". Such a wide 
range of spatial scales is not well sampled by any sin- 
gle VLA array configuration. The A array configuration 
at 1.4 GHz, for example, samples only spatial scales be- 
tween ~ 1.4' and ~l'/2. The usual approach of recovering 
the large scale structure by additional observations us- 
ing more compact array configurations is precluded by 
the fact that VLA configuration changes occur only ev- 
ery four months. 

I used a strategy for obtaining reliable, time-resolved 
images of an extended ob ject similar to that developed by 
IBietenholz et alJ l)2001alK Since a speed of c represents a 
proper motion of only ~l'/6 per month, any rapid evolu- 
tion must occur on the smaller spatial scales, which are 
well sampled by the A array configuration. The large 
scale structure as sampled by the B and C array con- 
figurations (spatial scales of 5" ~ 15'), therefore, is un- 
likely to change on timescales of ~ 2 months. I therefore 
made separate images of each A array session by supply- 
ing the maximum entropy deconvolution with a default 
image. The default image was made by imaging all the 
2003/2004 data, including that from both A, and the C 
and D array configurations, and subsequently smoothing 
the image by convolving with a Gaussian of FWHM 30". 
The use of the same default for both A-configuration im- 
ages will also serve to minimize any spurious differences 
between them since any differences of each from the com- 
mon default will be only such as are demanded by the 
data. 

3. RESULTS 
3.1. Images 

The full resolution image of 3C 58 is shown in Fig- 
ure n This image was made by combining all the data 
in Tabled including both A array configuration sessions. 
It therefore represents a time-average of the emission of 
the observing period of 2003 July to 2004 April. The 
position of the re centl y discovered pulsar, taken from 
iSlane et all l)2002|) and lCamilo et all l)2002|) . is marked. 
The image was made using robust weighting, with the 
AIPS robustness parameter set to 0, and the Gaussian 
convolving beam had a FWHM of 1'.'36. In addition, Fig- 
ure |5J shows a detai l of the central r e gion, including the 
"wisp" identified bv lFrail fc Moffetti l|1993j K I now turn 
to examining the variability of 3C 58's radio emission on 
several timescales. 

3.2. Expansion between 1984 an d £003 
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I examine first the longest timescales in order to de- 
termine the overall expansion rate. This can be most ac- 
curately determined from the pair of multi-configuration 
radio images from 1984 and 2003/2004 3 . I use the same 
approach to dete rmining the expansion as was used in 
iBietenholz et al.l 1)200 lbff . and repeat a brief description 
here for the convenience of the reader. The goal is to de- 
termine the overall or average expansion speed of the 
radio nebula. Since there are few well-defined, com- 
pact features, the expansion is measured not by deter- 
mining the proper motion of individual features, but by 
determining an overall scaling between a pair of images 
by least-squares. This was accomplished by using the 
MIRIAD task IMDIFF 4 which determines how to make 
one image most closely resemble another by calculating 
unbiased estimators for the scaling in size, e, the scal- 
ing and the offset in flux density, A and b respectively, 
and the offsets in RA and deck, x and y respectively, 
needed to make the second image most closely resemble 
the first. My chief interest is in the expansion factor, e, 
but because of uncertainties in flux calibration, absolute 
position, and image zero-point offsets caused by missing 
short spacings, all five parameters we re determined. This 
method was originally developed bv iTan fc Gull l|1985f) 
an d more details of i t s use in a similar situation are given 
bv IBietenholz et alJ l|1991bfl . Before running IMDIFF, 
the 2003/2004 image was convolved to the resolution of 
the 1984 one, namely 2"00 x 2"07 (FWHM) at p.a. -77°. 
The region exterior to 3C 38 was blanked so as to exclude 
it from the fitting. The fitting region encompassed 37.3 
square arcminutes or ^46,300 beam areas, and is indi- 
cated by a dotted line in Figure ^ The resulting expan- 
sion factor between the 1984 and the 2003/2004 images 
is 1.0027, with a formal uncertainty of ~ 0.0001. 

Although the formal uncertainty in e takes into account 
correlations of e with the other parameters (A, 6, x and 
y), it is still likely a significant underestimate because the 
true uncertainty will be dominated by systematic effects 
rather than by random noise. Systematic effects likely 
to influence my determination of e are problems in re- 
covering the largest scale structure in the deconvolution, 
bandwidth smearing, and the primary beam correction. 
I discuss each in turn. 

There are often problems in recovering the largest scale 
structure from intcrfcrometric observations such as these, 
and such problems might bias my determination of e. 
To eliminate this possibility, I high-pass filtered my im- 
ages using a Gaussian of FHWM 30", and obtained a 
marginally smaller expansion factor of 1.0022. To check 
whether biases are introduced by possible edge-effects, I 
determined e using only pixels above a certain cutoff in 
brightness. I determined e using the pixels with bright- 
ness > 3% and > 20% of the peak, resulting in fitting- 
areas of 90% and 52% of the above value, and values of 
e of 1.0021 and 1.0025 respectively. Similarly, tests show 

3 The image made from the 1991 data is not suitable for this pur- 
pose, firstly because it consists only of A-configuration data, hence 
does not reliably image structure larger than 1.4', and secondly 
because it suffers from bandwidt h smearing at distances > 45" 
from the phase-center (see e.g., ICottonl Tl999: Bridle & Schwab! 
I999j). This smearing can cause small but systematic radial shifts of 
the position of brightness peaks due to the non-symmetrical VLA 
bandpass. 

4 The IMDIFF program was slightly modified so as to treat 
blanking correctly. 



that the value of e is not sensitive to whether CLEAN or 
VTESS deconvolution is used. I conclude that my deter- 
mination of e is not sensitive to deconvolution errors at 
a level greater than 0.0006. 

I note that the total flux density of 3C 58 was found 
to be increas i ng; with time at a rate o f about 0.28% yr _1 
l)Greenlll987t lAller fc Revnoldslll985|) .. which would im- 
ply an increase of ~ 5.5% in the total flux density be- 
tween 1984 and 2003/2004. I find that the 5 GHz flux 
density of 3C 58 decreased by 0.7 ± 8.0% between 1984 
and 2003/2004 which is consistent with the earlier mea- 
sured secular increase but does not confirm it. 

The data of 1984 had a bandwidth of 12.5 MHz, and 
as a result will be subject to a small amount of band- 
width smearing at distances > 200" from the image cen- 
ter. Bandwidth smearing will smear the image radially. 
Such radial smearing might possibly bias the derived 
value of e, althoug h to first ord er it preserves surface 
brightness (see, e.g.. lCottonlll999ft . As a test, I averaged 
the 2003/2004 data over 25 MHz in frequency (from the 
original channel width of 6.25 MHz). Using IMDIFF to 
compare this artificially bandwidth-smeared image to the 
unsmeared original, I find the value of e to differ from 
unity by < 0.0003. As the amount of bandwidth smear- 
ing in this test was twice that in the 1984 data set, I 
conclude that any effect of bandwidth smearing of the 
1984 data set is small (and would cause my estimate of 
the expansion rate to be slightly but not significantly 
over-estimated) . 

Finally, different primary beam corrections might 
cause an error in the determination of e. I obtained my 
IMDIFF results from images uncorrected for the primary 
beam pattern of the VLA, which has the desirable result 
of making the noise in the image spatially uniform. If the 
two images were taken at a different observing frequency, 
the primary beam correction would differ, and a bias 
might be introduced. Since, however, the primary beam 
correction is not large, (being only ~ 7% near the eastern 
and western edges of the nebula), and since the 1984 and 
2003/2004 images are at almost the same frequency (the 
center frequencies differ by only 21 MHz), the difference 
in the primary beam correction between the 1984 and 
2003/2004 images is negligible. As a test, IMDIFF was 
run between images made separately from the two in- 
termediate frequencies in the 2003/2004 data set, which 
differ in frequency by a larger amount of 80 MHz. The 
resulting change in e was only 0.0005. 

As mentioned above, the uncertainty in e is likely dom- 
inated by difficult-to-determine systematic effects. The 
tests performed above, however, in no case resulted in 
a change > 0.0006 in the derived value of e, and I will 
adopt this value as a conservative uncertainty in e. 

The value of e implies 3C 58 's rate of expansion over 
the period 1984 to 2004 is 0.014 ± 0.003% yr" 1 . This 
value is consistent with, but more ac curate than, the 
one of .020 ±0.008% yr" 1 obtained bv lBietenholz et alJ 
( 2001b) from a variety of earlier, lower-resolution radio 
images. I note here that this expansion rate is much 
smaller than the rate of 0.124% yr -1 which would be 
expected for undecelerated expansion of 3C 58 since 
1181 A.D. 

3.3. Changes between 1991 and 2003/2004 
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Are there any changes in 3C 58's morphology 
apart from the relatively slow expansion? The 
Crab Nebula shows complex brightness changes, dis- 
tinct from the general e xpansion, which occur over a 
time-s cale of a few years llBietenholz et alJEM l2001at 
iBietenholz fc Kronberell992t) . Is there similar variability 
in 3C 58? The best pair images for addressing this ques- 
tion are those from 1991 and 2004, since they have the 
highest resolution and good signal-to-noise. In 1991, only 
A configuration observations are available, so only struc- 
ture < 1.4' in size can be imaged. This is not a serious 
limitation because any changes will likely relatively local- 
ized. A more serious limitation is that the bandwidth of 
43 MHz used in 1991 results in a considerable amount of 
bandwidth smearing — features far from the phase cen- 
ter will suffer radial smearing 45" of the phase center (see 
§ 13.21 above), which will limit the search for motions or 
changes to the region near the pulsar. The resolution of 
l'/36 and the time interval of 12.4 yr means I am sensitive 
to motions with projected speeds of > 1500 km s _1 . 

I formed a difference image by subtracting the 
2003/2004 image from that of 1991. To account for any 
possible small shifts in nominal position or differences in 
the amplitude scale etc., IMDIFF was first run to de- 
termine the best-fit flux density scaling and offset and 
RA and deck offsets (A, b, x, and y) between 1991 and 
2003/2004 images using only the portions within ±45" 
of the pulsar. The scaling in size (e) cannot be reliably 
determined from such a small region, so it was fixed at 
e = 1.0017, arrived at by assuming the expansion rate 
of 0.014 ±0.003% y^ 1 derived in §E3above. After ap- 
plying the relevant shifts and searings, I formed the dif- 
ference image. Because of both the IMDIFF scaling and 
the use of the same default image, this difference image 
should represent the smallest set of differences allowed 
by the data. 

This difference image is shown in Figure |3J There is a 
tantalizing hint of structure in the region near the pul- 
sar. The extrema within 20" of the pulsar are ±150 and 
— 120 /iJy bm _1 , or ±8% and —6% of the peak bright- 
ness. These extrema represent ~ 5<r, slightly larger than 
the extremes expected from a purely random variation 
over the whole image (which comprises n ~ 170,000 
beam areas). However, the distribution of image-plane 
noise is probably not strictly Gaussian, and the statis- 
tical significance of the extreme values in the difference 
image is therefore likely somewhat different than calcu- 
lated from Gaussian statistics. For this reason, the evi- 
dence for brightness variations near the center of 3C 58 
should be regarded as somewhat tentative. 

3.4. Variability over 110 Days 

Finally, I turn to determining whether there is any 
variability on the short time-scales of less than 1 year. 
For this purpose, I examine the differences between the 
two images made using each of the 2003 A configuration 
observations separately. In Figure H| I show a difference 
image formed by subtracting image from July 7 from that 
of August 9. No structure above the noise is visible in 
the difference image. Any changes in the radio morphol- 
ogy of 3C 58 over this period, therefore, were no larger 
than the noise. The rms brightness of the difference im- 
age over 3C 58 is 26 /iJy bm _1 , which is consistent with 
that expected from the background rms of the individ- 



ual images. The rms brightness over the region within 
±20" of the pulsar is not significantly larger. A limit of 
~ 160 //Jy bm _1 , or ~ 8% of the peak brightness can be 
set any changes over the period of 110 d 5 . 

In particular, the f eature identified as a "wisp" by 
iFrail fc Morfettl l|1993|) does not show any motion over 
this 110 d interval. I show profiles drawn in the R.A. 
direction through the wisp in Figure |5j To decrease 
the uncertainty, the profiles were averaged over 2" 4 in 
declination. Any changes between the two epochs are 
smaller than the profile uncertainty of ~ 15 ^tJy bm _1 . 
In particular, the wisp moves < 0"25 between these 
two epochs, implying a projected speed of < 0.05c or 
< 15, 000 km s- 1 (D = 3.2 kpc). 

4. DISCUSSION 

I produced new radio images of the supernova remnant 
3C 58. These new images are the highest quality radio 
images of 3C 58 so far, By comparing images made at 
intervals of ~20 yr, ~12 yr and 110 d, I was able to 
measure the expansion rate of the synchrotron nebula, 
and to establish upper limits on any variations or motions 
within the nebula. 

The new high-resolution images show a general ra- 
dio structure of 3C 58 very similar to that seen in ear- 
lier images at lower res olution and signal-to-noise (e.g., 
iRevnolds fc Alleil 11988) . Filamentary structure can be 
seen throughout the body of 3C 58, and the general tex- 
ture is very si milar to both tha t seen in the X-ray im- 
ages of 3C 58 (H ane et af]l2004fl and in radio images of 
the Crab Nebula (e.g., Bietenho lz et aUl2004j) . Several 
elongated loop-like structures are visible, particularly to 
the East. The narrowest filaments are likely to be un- 
resolved at my 1'.'36 resolution, but are relatively faint. 
The brighter filaments seem to have widths of a few arc- 
seconds. The most prominent filament is one near the 
pulsar, oriente d North-South, which was identified ear- 
lier as a wisp (Frail & Moffctt 1993). I will discuss it 
separately in § 14.41 below. 

4.1. The Expansion Rate and Age of 3C 58 

By comparing my observations to ones taken in 1984, I 
measured the current expansion rate of the radio nebula 
to be 0.014±0.003% yr -1 , corresponding to speeds 630± 
70 (Z)/3.2kpc) km s _1 along the major axis of the nebula 
and approximately half that along the minor axis 6 . This 
measurement of the expansion rate is consistent with, 
but considerably more a ccurate than the earlier one of 
IBietenholz et alJ l|2001b(l . As already pointed out in that 
reference, the low expansion speed of the radio nebula 
is corrobo rated by the slow proper motions of o ptical 
filaments (Fes en et ahl 119881: Ivan den Bergbl ll99Q). and 
by the low radial velociti es implied b y the low optical 
line-widths (<900 km s" 1 . lFesenlll983T) . 

5 As noted earlier, the computation of the exact statistical signifi- 
cance based on Gaussian statistics is likely to be in error. The given 
limits on the brightness variation over 110 d represent ~ 5 times the 
rms in the difference image. Tests performed b y adding random 
noise to the difference images of the Crab from Bietenholz et al. 
12004) show that the real brightness variations in the Crab are 
clearly distinguishable from the noise at this level. 

6 These speeds are calculated assuming symmetrical expansion 
about the center of the nebula. The observations constrain the 
expansion center to lie within about 1.5' from the pulsar position. 
The degree of one-sided expansion is not well determined. 
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All these measures of expansion speed of 3C 58 are 
much smaller than that expected if it were the undecel- 
erated remnant of SN 1181, in which case the expan- 
sion rate would be 0.124% yr -1 or ~ 5500 km s _1 along 
the major axis. Conversely, my measured expansion rate 
suggests an age for 3C 58 of ~ 7000 yrs with a 3<r lower 
limit of 4300 yrs, which is incompatible with the usual 
association with SN 1181. 

If 3C 58 is in fact the remnant of SN 1181, then, 
it must have be e n stron gly decelerated. As argued in 
iBietenholz et al.l l)2001b|) . to decelerate the synchrotron 
nebula, whose pressure is likely higher than that in the 
surrounding supernova ejecta, requires the ejecta to have 
been even more strongly decelerated. If the ejecta in- 
deed had been strongly decelerated, shell emission would 
be expected from their interaction with the interstellar 
medium. No shell emission, however, is seen in 3C 58 
in the radio despite several s earches (jBietenholz et all 
12001 Vi iRevnolds fc Alle3H§85|) . 

The simplest way out of this dilemma is to suggest 
that 3C 58 is not in fact the remnant of SN 1181, but 
rather is the remnant of a considerably older supernova. 
Several other lines of evidence suggest that 3C 58 might 
be rather older than 820 yr (many of these po ints have 
already been discussed in Ichevalierll2004 12005|) . 

1. The internal energy of the pulsar nebula, estimated by 
the synchrotron minimum energy, is larger tha n the cur- 
rent spindo wn luminosity multiplied by 820 yr l|Chevalierl 
1200412005]) . In order to have supplied the required neb- 
ular energy, the pulsar must already have spun down 
significantly, which is unlikely for a pulsar only ^820 yr 
old with a spindown age of 5380 yr. 

2. An age of ^820 yr implies rapid expansion of the PWN 
since the supernova, which requires an impla usibly low 
densi t y in the surrounding supernova ejecta (Chevalier 
12004 |2005|) . The expansion of the PWN into the su- 
pernova ejecta is driven by the spindown energy of the 
pulsar, which in the case of 3C 58 is ^20 times lower 
than that of the Crab nebula. Assuming similar super- 
nova ejecta in both remnants, one would expect 3C 58 to 
expand considerably more slowly than the Crab nebula, 
yet if 3C 58 is only ~ 820 yr old, then its average expan- 
sion speed is roughly double that in the Crab nebula. 

In addition, the mass swept up by the expanding PWN, 
estimated from the thermal componen t of the X-ray 
emission near the edge of the nebula (Boc chino et alJ 
1200 ll ISlane et all |2004) , is > 30 times larger than ex- 
pected for an age of ~ 820 yr and a spindow n luminosity 
of 2.7 x 10 37 erg s" 1 l)Chevalierll2004 12001 . 

3. The temperature of the thermal X-ray emission sug- 
gests a shock velocity of ^340 km s _1 , whereas if the 
PWN is expanding into the freely expanding supernova 
ejecta with an age of only ~ 820 yrs, a shock velo city of 
~1000 km s- 1 is expected l|Chevalierll200l ElO^ . 

4. The spindown age of the pulsar , P/2P , is 5380 yr 
ijMurrav et alJ l2002t iRansom etafl l2004|) . Although 
spindown age is n ot a reliable measure of true age (e.g., 
iKasoi et al.ll200TI) . the spindown age of the pulsar is 
nonetheless consistent with an age of several thousand 
years for 3C 58. 

5. The very sharp and rather low frequency spec- 
tral break observed in 3C 58 ijWoltier et all Il997t 



iGreen fe Scheuerl H992) suggests an age older than ~ 
820 years. 

I note also that the X-ray observations show that the 
neutron star is much cooler than expected for a normal 
cooling curve and an age of ~ 820 yr, although vari- 
ous "non-standard" cooling mechanisms or a high-mass 
neutron star could likely explain this discr epancy (e.g., 
IPage et al.ll200l lYakovlev fc Pethickll200l . 

4.2. Alternative Remnants of SN 1181 

If 3C 58 is not the remnant of SN 1181, an obvious 
question arises: where is the remnant of SN 1181? At 
only ~820 years old, it might be expected to be fairly 
bright and have an angular size of > 1'. Why has it not 
been discovered? The Chinese, Korean and Japanese 
records place SN 1181 within ~1° of the Galactic- 
coordinates I = 130°, b = +3° iSteohensonl Il971i 
iClark fc Stephenson! Il977l iStephenson fc Greenl 11999^ 
All the other historical supernovae were at distances 

< 4 kpc, and since SN 1181 was not unusually faint, 
it was likely within 5 kpc. 

The most sensitive radio surveys covering the re- 
gion are the Cana dian Galactic Plane Survey (CGPS; 
iTavlor et al] l200l and th e NRAO VLA Sky Survey 
(NVSS: lCondon et alJ lf998). No obvious candidate with 
a size > 1' is visible on either survey, with the 3cr limits in 
surface brightness at 1 GHz, Sighz, being ~ 15xl0~ 22 
and ~ 8 x 10~ 22 W m~ 2 Hz -1 sr" 1 , respectively 7 . In 
addition, the sensitive targeted searches for 3C 58 's shell 
emission would be expected to have shown any emission 
fro m other remnants nearby: iRevnolds fc Allerl l)1985fl 
and IBietenholz et alJ l)2001b|) placed limits on Sighz 
of ~ 5 x 10~ 22 W m -2 Hz -1 sr" 1 . In particular , the 
327 MHz observations of IBietenholz et all l|2001b|) had 
a field of view of ~1.2° and thus would be expected to 
have detected the remnant of SN 1181 if it were suffi- 
ciently bright. 

If 3C 58 is not the remnant of SN 1181, we can there- 
fore say that the latter's remnant must have a radio sur- 
face brightness < 5 x 10~ 22 W m~ 2 Hz -1 sr _1 . This 
surface brightness is an order of magnitude fainter than 
the faintest known shell remnants. Is it reasonable to 
suppose that the remnant of SN 1181 might be so faint? 

There are two reasons to think that such a faint rem- 
nant is not improbable. Firstly, estimates of the Galactic 
supernova rate suggest that ther e were 12 to 24 Galactic 
super novae in the last 1000 yr ijReedl 120051 iCappellarol 
l2003HStrolgerll2003|) . One would therefore expect to ob- 
serve the same number of supernova remnants with ages 

< 1000 yr, which in the radio are likely to be visible any- 
where in the Galaxy. Currently, despite several searches, 
only about 10 know n remnants have ages or estimated 
ages < 1000 yr (e.g . . lGreenl | 2004t iMisanovic et~alll200l 
iSramek et alJ Il992t IGreen! 1 19851) . The implication is 
that there are likely ~ 2 to 14 as yet undetected young 
remnants. It is therefore plausible that the remnant of 
SN 1181 is one of these as yet undetected remnants. 

7 To facilitate inter-comparison, I cite surface brightnesses at 
1 GHz, and convert measurements at different frequencies using 
an assumed spectral index of —0.5. For example, both the CGPS 
and the NVSS surveys were at 1.4 GHz. Note also that some 
artifacts due to the presence of 3C 58 are visible in both the CGPS 
and NVSS surveys, hence the present surface brightness limits are 
somewhat higher than the nominal values for those two surveys. 
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The second reason is to think that we may not yet have 
discovered the remnant of SN 1181 involves a supernova 
of almost the same age: SN 1054, which is one of only six 
known historical supernovae. The only presently known 
radio emission associated with SN 1054 is the Crab neb- 
ula and pulsar. The nebular radio emission can be traced 
entirely to the Crab pulsar, which has a very high spin- 
down luminosity. No radio emission from the supernova 
shell has yet b een detecte d : the most sensitive search to 
date being by iFrail et all ^)!)5i , which found no radio 
emission from the supernova shell with a 3<r upper limit 
on Sighz of < 5 x 10~ 22 W m~ 2 Hz^ 1 sr^ 1 . This upper 
limit is fainter than the detection limits of the CGPS 
and NVSS surveys mentioned earlier. Therefore, if it 
were not for its energetic pulsar, the radio remnant of 
SN 1054 might not yet have been discovered. It is then 
not unreasonable to suppose that SN 1181 produced a 
similarly faint shell, but did not produce an energetic 
pulsar. 

A possible, but unlikely, scenario would be that 3C 58 
is the composite remnant of a binary where both stars 
have undergone supernova explosions, with the earlier 
supernova leaving the slowly expanding remnant, and the 
later, low energy one providing the historical supernova 
event in 1181 A.D. This scenario seems rather contrived, 
since having two supernovae within a few thousand years 
of one another is quite unlikely 

4.3. Rapid Variability near the Pulsar 

No conclusive evidence is seen for rapid variations, ei- 
ther motions or changes in the radio brightness, near 
the pulsar. The upper limits on rapid variability are 
~ 8% of the peak brightness over a period of 110 d. In 
contrast, in the case of the Crab Nebula, a distinct pat- 
tern of elliptical ripples is seen in radio difference images 
over similar time intervals. However the relative ampli- 
tudes of the difference r ipples in the Crab is of only 6% 
ijBietenholz et al J 12004 l2001ajh slightly lower than the 
upper limit on the variations occurring in 3C 58. 

4.4. The Wisp 

As mentioned, the feature that IFrail fc Moffettl (1993) 
identified as a "wisp" is clearly seen (see Figs, an d 
EJ), but doe s not appear to move over 110 d (projected 
v < 0.05c). IFrail fc Moffettl l|1993D interpreted the wisp 
as being associated with the pulsar wind termination 
shock. Since 3C 58's wisp, unlike those in the Crab neb- 
ula, does not appear to move rapidly, the question arises 
of whether it is in fact associated with the pulsar wind 
termination shock. 

The wisp's morphology does not argue strongly for an 
association with the wind termination shock, as filaments 
just as narrow and in a variety of orientations are seen 
in the remainder of the nebula. Unlike the wisps in the 
Crab, 3C 58's wisp is relatively straight and does not 
unambiguously suggest a circular feature seen in projec- 
tion. 

The spindown luminosity of the pulsar is now known 
to be ~ 2 7 x 1 37 erg s" 1 ijMurrav et alJ 120021: 
iRansom et all 12004 . therefore the approximate loca- 
tion of the termination shock can be calculated by bal- 
ancing the ram pressure of the pulsar wind (EAirr 2 ) 
against the nebular pressure. Standard synchrotron 
arguments give a minimum pressure in the nebula of 



4 x 10~ 10 dyne cm -2 , so the termination shock radius 
should be at a distance of < 4 x 10 17 cm, corresponding 
to <8" for D = 3.2 kpc, with the distance being lower if 
equipartition does not obtain and the pressure is higher 
than the synchtrotron minimum pressure. 

The wisp is at an angular distance of only 3'.'5 from the 
pulsar, somewhat closer than the expected termination 
shock, although, as mentioned above, this might be due 
to departures from equipartition. Moreover, if the wisp 
is in fact a circular feature, then its inclination angle is 
> 60° and its radius corresponds to > 7", compatible 
with an association with a termination shock. I note 
that, in this case, the space velocity of the wisp would 
be considerably higher than the measured projected ve- 
locity, perhaps even as high as the velocities seen in the 
Crab wisps. 

High-resolution X-ray images of 3C 58 in fact sug gest a 
highly inclined torus with a jet (|Slane et alJ l20041. The 
radio wisp is roughly coincident with the western edge 
of the X-ra y torus (as has already been pointed out by 
lSlandl2005|) . although it is possibly slightly larger. 

The location, orientation and relative prominence of 
the radio wisp and its similarity to the X-ray "torus" , 
therefore, argue strongly that 3C 58's radio wisp is indeed 
associated with the wind termination shock despite its 
lack of apparent motion. 

5. CONCLUSIONS 

New VLA radio observations of 3C 58, and reanalysis 
of older observations from 1984 and 1991 indicate that: 

1. A filamentary structure is present throughout the 
body of the nebula. 

2. The radio nebula expands at a rate of 0.014 ± 
0.003% yr" 1 . This is much slower than expected if 3C 58 
is the undecelerated remnant of supernova 1181 A.D, but 
suggests instead an age of several thousand years. 

3. The fractional variability in brightness (other than 
that due to expansion) is < 8% both over a period of 
110 d and over the period 1991 to 2003/2004. There 
is inconclusive evidence for brightness variations of a few 
percent in the central portion of the nebula over the latter 
period of approximately a decade. 

4. A linear feature, called "wisp" very near the pulsar, 
is likely related to the pulsar outflow. 

5. The wisp does not appear to move over a period of 
110 days, its projected speed being < 0.05 c. This stands 
in contrast to the Crab nebula, where rapid motion of 
the wisps with speeds of ~ 0.3 c are seen. 
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TABLE 1 

New Observing Sessions in 2003 - 2004 a 



Date 


Array Configuration 


Observing Time b 




(hrs) 


2003 Jul 7 


A 


12 


2003 Aug 9 


A 


12 


2003 Dec 30 


B 


4.5 


2004 Apr 19 


C 


3.5 



a All observations made in spectral line mode, using two 
intermediate frequencies (IF), centered at sky frequencies 
of 1.4649 and 1.3851 GHz, each with 7 spectral channels 
of width 6.25 MHz. 

b The total observing time, including calibrator observa- 
tions 
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Fig. 1.— An image of 3C 58 at 1.4 GHz, made using all the A, B and C array-configuration VLA data from 2003 and 2004. The FWHM 
of the Gaussi an convolving beam was l'/36. The c ross marks the position of the pulsar, PSR J0205+6449 (RA = 02 h 05 m 37.92!, deck = 
64° 49' 42 / /8: ISlan"e et alj|200l ICamilo et alj|2002l) . which is known to about 0'.'5. The grey-scale is labelled in mjy bm _1 , and the rms 
of the background was ~ 12 fj,]j bm" 1 . The background sources, including the bright extragalactic double at R.A. = 2 h 6 m 31?7, deck = 
64° 54' 14 " 5, w ere subtracted from the visibility data before the final maximum-entropy deconvolution, hence are not visible in this image 
(see text S I3.1l for details). The dashed line indicates the fitting region used for the expansion calculation (see § 13.21 . For a detail of the 
region near the pulsar, see Figure |21 
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Fig. 2. — An image of the region near the pulsar in 3C 58. This figure shows the central part of the image in Figure The cross again 
marks the position of the pulsar, known to ~ 0.5". 
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Fig. 3. — D ifference between the images from 2003/2004 and 1991, after the best-fit expansion, translation and scaling in brightness 
(see text S 13.31 , Due to the bandwidth smearing in 1991 image, the difference image becomes less reliable beyond ~ 45" from the pulsar 
position, which is marked with a cross. The FWHM of the Gaussian convolving beam was l'/36. The contour, included for reference, is 
the 30% contour of the 2003/2004 image in Figure IT1 The grayscale is labelled in fijy bin -1 . 
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Fig. 4. — Central portion of the difference image between 2003 July 7 and August 9, i.e. between images made using separately each 
of the 2003 A configuration data sets. The FWHM of the Gaussian convolving beam was 1'.'36. The rms variation in the difference image 
is 25 fj,3j bm" 1 . The plotted contours are the 2 and 20% contours from the total intensity image in Figure shown for reference. The 
position of the pulsar is marked by a cross. 
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Fig. 5. — Profiles in R.A. through the "wisp" on 2003 July 7 and 9. The profiles are at a declination of 64° 49' 40"2 (J2000). To improve 
the signal-to-noise ratio, the images were smoothed in the declination direction by boxcar averaging over 2'.'4. The FWHM Gaussian 
resolution in the R.A. direction is 1'.'36. The uncertainty is ~ 0.015 mjy bm" 1 in each profile. 



